MONOTONE UNITARY FAMILIES 

^. . DANIEL GRIESER 

O 

O , 

^sj Abstract. A unitary family is a family of unitary operators U(x) acting on 

a finite dimensional hermitian vector space, depending analytically on a real 
r parameter x. It is monotone if iU'(x)U(x)~ 1 is a positive operator for each x. 

We prove a number of results generalizing standard theorems on the spectral 
S { theory of a single unitary operator Uo, which correspond to the 'commutative' 

case U(x) = e tx Uo- Also, for a two-parameter unitary family - for which there 
Ui , is no analytic perturbation theory - we prove an implicit function type theorem 

for the spectral data under the assumption that the family is monotone in one 

argument. 

< 

^ ' 1. Introduction 

a'. 

Let U(x) be a family of unitary operators on a Hermitian vector space V of 
dimension M < oo, depending real analytically on x € R (or an interval in K). 
Then 

(1) D(x) := -U'^Uix)- 1 

VO . is symmetric. Here U'(x) is the derivative with respect to x. We call U monotone 

if D(x) is a positive operator for all x. Denote 

W(x) = Ker(I - U{x)) and Z = {x : W{x) ^ {0}}. 

Thus x G Z iff U(x) has eigenvalue one. 

A model case for this setup is U(x) — e zx Uo for a unitary Uq. Then W(x) is the 
eigenspace of Uo with eigenvalue e~' lx . Standard facts from the spectral theory of 
Uq may be restated in terms of Z and W(x), for example: 

• Z is a 27r-periodic sequence, having exactly M terms in each half-open 
H | interval of length 2ir (counting 'multiplicities'). 

• If I is an interval of length less than 2ir then the spaces W(x), x £ I, are 
linearly independent (even pairwise orthogonal). 

• If ||(7 - U(x ))ip\\ < s\\ip\\ for some <p £ V \ {0}, x e K and e > then 
e~ lx ° lies within distance e of an eigenvalue of Uo (and so dist(xo,-Z) < 
ne/2). Furthermore, if e' > e and P denotes the orthogonal projection to 

® X :\e-**-e-**o\ <e 'W(x)then 

(2) l|p--fVII<plM|. 

In this paper we prove generalizations of these facts to arbitrary monotone uni- 
tary families, see Theorem Q] in Section [5] and Theorems [5] and in Section [5] 
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The estimates are expressed in terms of uniform bounds on the first and second 
derivatives of U: Assume d m in, ^max, ^2 > are such that 

(3) d min I < D(x) < d max 7, ||J7"(s)|| < d 2 for all x. 

Of course, such constants exist always locally. However, in applications of our 
results, see [2], it is essential that the estimates are uniform in terms of ((3]), and do 
not depend on additional data like separation of elements of Z, see the explanation 
below. 

An important difference between a general unitary family and the model case is 
that U(x) and U(x') do not in general commute for x ^ x' . A consequence of this 
is that, while one may take a logarithm of U, i.e. find an analytic family A{x) of 
symmetric operators such that 

(4) U(x) = e lA{x) for all x, 

it is usually not true that D(x) equals A'{x), or even that positivity of D(x) implies 
positivity of A'(x). The opposite implication is true, however. See Section [U 

Finally, we prove a result on two-parameter perturbation theory. Recall the 
main result of one-parameter perturbation theory (see [4], [3]): For an analytic 
family U(x) of unitary operators on V, there are real analytic functions Hj(x), 
<fij{x) for j = 1, . . . ,M = dimy having values in R and V, respectively, such that 
for each x the eigenvalues of U(x) are e l,J " j ( x \ with corresponding orthonormal basis 
of eigenvectors ipj (x) : 

(5) U(x)(p j (x) = e i,i *W(p j (x). 

It is well-known that the analogous statement for two-parameter families of opera- 
tors is false in general. However, we prove a related implicit function type theorem 
for the spectral data of a two-parameter family for which the dependence on one 
parameter is monotone. It may be regarded as the natural unitary family gen- 
eralization of the one-parameter perturbation theory. See Theorem [7] in Section 

El 

The analytic functions \ij and (fj play a central role in the proofs of our theorems. 
It is essential to control their derivatives. For /ij this is easy from ([3]). However, 
ipj may vary wildly whenever e t>lj is very close to another eigenvalue. To control 
this variation is the main technical problem in the proof of the generalization of 
(f2|), Theorem [3] Note that, unlike in the model case, one may not assume that 
the e Vi ( x ) for fixed x but varying j are uniformly separated, or equivalently that 
the elements of Z are uniformly separated. This can already be seen in the simple 
example U(x) = e lxL , where L is a diagonal matrix with positive diagonal entries 
that are independent over Q: If M > 2 then for any e > there are x, x' G Z 
satisfying < \x — x'\ < e. 

The problems we study here arose in the context of a singular perturbation 
problem: In [2] we study the eigenvalues and eigenfunctions of the Laplacian on a 
space X N which has a fixed compact part connected by cylindrical necks of length 
N > 0, and in particular their asymptotic behavior as TV — > oo. The unitary families 
arise from the scattering matrix of the limit problem (infinitely long 'necks'). 
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2. Eigenvalue distribution 



Theorem 1. Let U be a monotone unitary family on K. Then Z cl is a discrete 
subset, and more precisely for all A < B 



(G) 



^2 dimW(x) - — / trD(x)i 

A ^ ~ ^ D J A 



B 

/ dx 

. ^ I A 

x:A<x<B J A 



<M(:=dimV) 



In the special case U(x) — e zxL Uo, L > 0, this implies the asymptotics 

trL 

~2tT 



Y^ dimW(x) ~—B, B 

x:0<x<B 



This is the Weyl asympotitcs of a quantum graph, see for example [T] . We give a 
much simpler proof than [T] . 

First, we differentiate |(5J) and obtain monotonicity of the functions fij from 
monotonicity of U: 

(7) $ = (ilTU-^Vj) = (Dvj,^) > 0. 

I 

Also, for each x we have 

(8) W(x) = span{ipj(x) : (ij(x) £ 2nZ} 
Proof of Theorem [7J We have 

J2(Hj(B)-Vj(A)) = / 22(D(x)cp j (x),ip j (x))dx 

J A 



3 
B 

tr D(x) dx 

I A 

Since the jj,j are strictly increasing, we get from © 



fi 3 (B)-^(A) 



J2 dimW(x) =Y#{keZ: N {A) < 2nk < Mi (*)} = £ f W£j_JWV + R . 

x:A<x<B j j 

with \Rj\ < 1, and this gives ©. D 

3. ElGENSPACES 

In this section we consider monotone unitary families satisfying the estimates 
© . By J7]) we have 

(9) d mia < fi'^x) < d max for all x and j = 1, . . . , M. 

First, we have independence of cigenspaces. 

Theorem 2. Let U be a monotone unitary family satisfying ©. Let L be an 
interval of length at most 2 f"l"r ■ Then the spaces W(x), x £ /, are independent, i.e. 

(10) If ip x £ W{x) for each x £ / n Z and 2_, 'fix = i/ien y^ = Vx. 
The following theorem gives a stable version of almost orthogonality. 
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Theorem 3. Let U be a monotone unitary family satisfying ([3|). Assume ^el^\0 
satisfies 

(ii) \\(i-u(x )M<4vl 

Then 

(12) d ist(zo,£)<^-7^. 

Furthermore, there is a constant C only depending on d m i„, d max , d2, M such that 
the following holds: Suppose < e < e' < C _1 and e/e' < C _1 . Denote by Pw the 
orthogonal projection to W — @\ x _ x i< e > W(x). Then 

(is) y-p w v\\<c^)^M\. 

In the proofs of Theorems [2] and [3] we will need the following estimate, which 
replaces orthogonality of the eigenspaces of a unitary operator. 

Lemma 4. If if € W(x), ip S W{y) and x ^ y then 

(14) \(D(x)^)\<^\x-y\-M-\m. 

Proof. By Taylor's formula, 

U{y)U{x)- x =I + (y- x)U'{x)U{xr x + (y - x) 2 R, \\R\\ < ^, 
so if U(x)<p — ip, U(y)tp — if) then 

foV) = (J7(s)-V,^(v) _ V> = (UivMx)- 1 ^) 

= (</>,V>) +i{y-x){D{x)ip^:) + (y-x) 2 (Rip,i{;), 
and this gives (fT4")l . D 

Proof of Theorem^ Let ^ x G W(a;) for ir e / fl Z, and assume ^ x <^ x = 0. Let 
ip Xo have maximal norm among all <p x . Then = (D(xo)<p Xo , ^2 X <p x ) gives with ([3]) 
and IE 



rfminll^oll 2 < ( D {xo)<Px ,<Px ) = 



x^x 



< s; *i j i • nv«oii • iiv-ii < (Af-i)*i/i • n^oii 2 , 

X^XQ 

so if d m i n > (M — l)4f |I| then <p Xo = and hence ip x = for all x. This implies 
the claim. □ 

Proof of Theorem^ The first estimate follows easily from the fact that, by the 
lower bound in @, an eigenvalue close to one of U(xq) will turn into an eigenvalue 
equal to one of U(x), for some x close to xq\ Let B{x) = I — U{x) and let \j(x) — 
l_ e Vj(z) be the eigenvalues of i?(x). The assumption (fTTj) implies that |Aj(xo)| < £ 
for some j, and this implies dist(/ij(xo)j 27rZ) < §£, and then //■ > d m i n shows that 
there is an x satisfying \x — xq\ < irs/2d m in and [ij(x) € 27rZ, hence i e Z, so (fT2|) 
follows. 
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For 5 > let Ps(x) denote the orthogonal projection to the sum of the eigenspaces 
of B{x) with eigenvalues |Aj(x)| < 5. Then ||i?(xoV| < e\\<p\\ implies 

(15) \\<p - Ps(x )cp\\ < 6 -\\<p\\ 

(see (|5J), which also applies to normal operators). To make this a good estimate, 
we want to take 5 >> e. Our goal is to replace Ps{xq) by Pw here. The idea is 
that eigenspaces of B(xq) with eigenvalue |Aj(xo)| < 6 will turn into nullspaces of 
B(x) for some x within 25/d min of Xo, by the first part of this proof. However, the 
variation of eigenspaces is much less well behaved than the variation of eigenvalues: 
An eigenspace may change rapidly with x if the eigenvalue is very close to another 
eigenvalue. Therefore, we need to consider not single eigenspaces but rather clusters 
of eigenspaces. 

The variation of eigenspaces is given as follows (see [3]): Fix x. If B(x) has no 
eigenvalue on the circle |A| = 5 then, with a prime denoting derivative in x, 

(!6) P's= E E Y^YiPjB'Pu + P.B'P,). 

j:|Aj|<5 fc:|A fc |>5 J k 

Here, all quantities are evaluated at x, and Pj is the orthogonal projection to 
spant/3j. Taking norms and using orthogonality of the Pj one obtains from this, 
using | IB' 1 1 < d max , 



(17) \\Ps\\ < d max M ( min [A,- — Afc 1 

\|Aj|<5, |A fc |><5 

We need to choose 5 carefully to make the spectral gap not too small: Let s — 
(e' /e) 1 '( m+1 > and consider the M disjoint subintervals 

[es k : es k+1 ) ioik = l,...,M 

of (e,e')- Since B(xq) has M eigenvalues and one of them has absolute value < e, 
at least one of these intervals contains no |Aj(xo)|. Assume 

(18) [8,Ss), S — es k , contains no |Aj(xo)|. 

We then have: 

A) The eigenvalues of B(xq) with |Aj(xo)| < S are in 1-1 correspondence with those 
x G Z (counted with multiplicity dimiy(x)) satisfying 

I \<rx> 2S 

\x -xo\ < o := - — . 

(Proof: Each such eigenvalue turns into a zero of Aj(x) for such an x, by 
the argument at the beginning of this proof. Conversely, if Aj(x) = then 
|Aj(x )| < |x-xo|dmax since |A^| = fij < d max and hence |Aj(xo)| < 2^^^- < 6s 
provided e/s' is sufficiently small (and therefore s big), and by (fT8|) this implies 
further |Aj(xo)| < 5.) 

B) The smaller interval (S + 2<5^ mai , 5s — 25 ^^ ) contains no | Aj (x) | for any x with 
|x — xo| < 5' . (Follows directly from |A'-| < rf mia -) 

The length of the interval in B) is 5{s — 1 - 4^^) which is > 5s/2 > if s 
is sufficiently big. Choose 5 in this interval, then we get from (|17[) ||Pj(x)|| < 
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d m axM {8s/2)~ for \x — xo\ < 5'. Integration gives 

(19) \\P-( x )-p.^ )\\< £l ~t^ M iov\x-xo\<S'. 

This implies 

(20) \\1> - Pt(xo)1>\\ < ei\\1>\\ iorx£Z, \x-x \<S', ip£\V(x) 

since Pg(x)tp = tp then and Pg(xo) — Ps(xq). 

Next we want to extend this estimate to ip £ W :— @\ x _ x i<j/ W(x). For this 
it is essential that, by (fT4|) and the positive definiteness of D(x), the angles between 
different W{x) are bounded away from zero. To carry this out, we first derive from 
(fT4|) an estimate where all D(x) are replaced by D(xq)'- From D' = f/"[/ _1 + D 2 
we have \\D'\\ < d 2 + d max 2 ; integration yields ||D(x )|| < ||-D(x)|| + (d 2 + d max 2 )S' 
for \x — xq\ < S' , and then (TTJJ gives, with Dq := D(xq), 

(21) |(^o^,W>l<*"ll^*H-|WI. 5" :=5'{d 2 +d max 2 +d 2 ) 

for ij) x £ W(x), ip y £ W(y), x ^ y and \x — Xo\ < S', \y — Xq\ < §' . Introduce the 
scalar product (ip,ip)r) := (Dotp,ip) on V, with norm HV'H.Do = y^flo^Pi W)> then 

(22) dminUV'll 2 < IIV-IIa, < dmaxllV-f 

so ([2Tj) gives 

(23) \(^,%)d \ < 1 — UAd ■ Uy\\ Do 

for the same tjj x ,ijjy as there. By simple standard calculations this implies 

(24) || ]T i> x \\l >(l- -f-(M i)) £ IIiMa, 

\ flmin / z 

where the sums are over all x £ Z with |aj — Xq\ < 6', tp x £ W(x) are arbitrary 
and M is the number of summands. Now by A) above M < M . Since 6 < e' the 
expression in parantheses is > | for sufficiently small e' , so (|24| gives J^ ||^ x [[x) < 
2|| E x V>*IId , which with m gives £ x ll^ll 2 < 2^|| E.^ll 2 and so 



r- ^max 



(25) ^||^||<W2A/-p||^^l|. 

V &min 

x x 

We return to |20|). If V £ W = (& lx _ Xol < s > W(x), ip = J2 X i>x then we get from 

(USD 

||V - P S {x )ip\\ < J2 HVfc - Ps(xo)ip x \\ <£iJ2 H^*ll 



(26) <£2||V||, e 2 :=e 1 J2M-^ 

This means that W is close to Wo : = RanP{(xo). Now by A) above, W and Wo 
have the same dimension, so this implies by standard arguments that Wo is also 
close to W, more precisely, with Py/i : V — > W the orthogonal projection, 

(27) U - PwH < e 2 U\\, ^£W . 
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Finally, assume that ip E V satisfies (fTTj) . Then (fl"5|) with 5 as above together with 
(|27|) give, again by standard facts about distances of subspaces, 



(28) ||^-/WII<(£2 + |)IM|. 

Putting everything together, we have S — es k for some k E {1, . . . , M} with s — 
(e'/e) 1 ^ M+1 \ which implies e/S < s' 1 and 6' < e' . Also, e x = ~|^f and e 2 = 

£ i W / 2Af|^. Altogether, the right hand side of (28]) is bounded by Cs _1 [[^|[, 
and since the left hand side only decreases when replacing W by the bigger space 
®\x-x \<e' W(x), the Theorem is proven. □ 

4. MONOTONICITY OF U AND OF ITS LOGARITHM 

Denote by S(V) and U{V) the spaces of symmetric resp. unitary operators on 
V. The map S(V) — > U(V), A *— » e lA has non-singular differential everywhere and 
is surjective, so it is a covering map. Hence any curve U : x i— » t/(x) in U(V) may 
be lifted to a curve a; i— ► A(x) in 5(F) (that is, £/(x) = e lA ^ for all x) and the 
lift is unique if one prescribes it for one value of x. Furthermore, the lifted curve is 
analytic if U is. 

Proposition 5. Let x t— » A(x) be a C 1 family of symmetric operators and U(x) — 
e lA( - x \ Then 

1 f 1 

(29) -U'U- 1 = e lTA A'e- lTA d T . 

i Jo 

Here, a prime denotes differentiation with respect to x, and U,U' , A, A' are taken 

at a fixed x. 

Proof. Let W(t,x) = j^e itA{x \ Then 

® W = -^-^-e UA = ^-Ae UA = A'e ltA + iAW. 
at i ox at ox 

\l a itA 



Now fix x, and let B(t) = A'e 1 . The solution of the ordinary differential equation 

di ' 



4iY(t) = B{t) + iAY(t) with F(0) = is 



V(t)= / e l{t ~ T)A B{T)dT 
Jo 

as can be verified directly. (This is called Duhamel's principle.) Now W(-,x) 

satisfies the same first order ODE and initial condition as Y, so it follows that 

W(t, x) = Y(t) for all t. Now rewrite Y(t) = / Q * e lrA B(t - r) dr and set t = 1 to 

obtain (HH). □ 



Corollary 6. Let x i— ► A(x) &e a C 1 family of symmetric operators. If A'(x) > 
/or eac/i x t/ien i/ie unitary family U(x) — e ^ x > is monotone. 

Proof. Positivity of A' implies positivity of e lTA A'e~ lTA for each r, so the claim 

follows from dH). □ 



The converse is not true. As an example let A = I j , B = I 

. , „ , - , a , \ * „ ^, .-^ 4 /cos 7tt — sin 7tt\ . 

with < b < 1 and A(x) = An + xB. Then e 2rj4 ° = . is rotation 

V ' V SUl TTT COS 7TT / 
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by itt, and a short calculation shows that L e %TA °Be~ %TA dr = ^-^-1 (this is also 
clear without calculation since the result must be rotation invariant with trace 
equal to tr B = 1 — 6; in essence, the negative direction of B gets averaged away 
against the positive direction). Therefore, U(x) = e tA ^ is monotone near x = 
but A'(Q) = B is not positive. 

5. TWO PARAMETER FAMILIES 

Theorem 7. Let U(x,y) be a unitary operator in a finite- dimensional Hermitian 
vector space depending real analytically on x,y € R. Assume 

(30) i|^7-i > at(x ,y ). 

Then the set {(x,y) : U{x,y) has eigenvalue one} is, in a neighborhood of (xo,yo), 
a union of real analytic curves x = Xj(y). The corresponding projections Pj(y) 
to the eigenspace of U(xj(y),y) with eigenvalue one are also analytic functions 
°f V ¥" Ho > extending analytically to y = yo, and J2j Pj( y o) * s the projection to 
ker(J- U(xo,yo)). 

Note that in general it is not true that the eigenvalues and eigenprojections of 
U(x,y) may be arranged as real analytic functions of (x,y), see (3 , II. 6.1. While 
the example given there (in the analogous case of self-adjoint operators) does not 
satisfy the positivity assumption (|30|) . it can be easily modified so it does, by 

adding a multiple of the identity. Explicitly, one may take A(x, y) = I I and 

U(x,y) = e iA ^ and {x 0l y ) = (0,0). 

Note also that the statement of the theorem reduces to the well-known facts 
of one-parameter perturbation theory in case U(x,y) = e lx U(y), for an analytic 
one-parameter family of unitary operators U(y). 

Proof. Let w.l.o.g. x = yo = 0. 

We first consider the case [7(0, 0) = /. Let A — klogU near (x,y) = (0,0). 
Then the operators A(x, y) are self-adjoint, A(0, 0) = 0, and dA/dx(0, 0) > since 
it equals j^U -1 ^^) by ([29]) . and we need to prove that the set S = {(x,y) : 
A(x, y) is not invertible } is a union of real analytic curves as claimed \j 

If A(x, y) = xA + yB is linear in x, y, then (since A > 0) A and B may be 
diagonalized simultaneously, hence may be assumed to be diagonal, and then it is 
obvious that S is a union of lines, Xj(y) = ybj/aj, where aj,bj are the diagonal 
entries of A, B, respectively. In general, write A(x,y) = xA + yB + C{x,y) with 
C(x,y) — 0(\x,y\ 2 ) and w.l.o.g. A, B diagonal. Then, if the dimension of the vector 
space is M, 

M 

det A(x,y) = Y[(xa 3+ yb J ) + 0(\x,y\ M+1 ), 

3 = 1 

and a standard argument (using polar coordinates) shows that the zero set of this 
function is a union of real analytic lines x — Xj(y), having tangents xa,j + ybj = 
at the origin. 

If U(Q, 0) is arbitrary, let W = U — I (where / denotes the identity) and Vo = 
Ker 11^(0, 0) and V\ its orthogonal complement. Let Wki(x,y), k,l — 0,1, be the 



I am grateful to Y. Colin-de-Verdiere for a fruitful discussion on this 
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'submatrices' of W(x,y) corresponding to the decomposition Vq © V\. Then Woo, 
Woi and Wio vanish at (x,y) ~ (0,0), and W\\ is invertible at (0,0) and hence in 
a neighborhood. Then the equation W(vq © v\) = 0, where vo <E Vo, vi G Vi, is 
equivalent to W'v = 0, where W = Woo — WoiW^Wio, and «i = -W^W w v . 
Therefore, U{x, y) has eigenvalue one iff the operator U'{x, y) — W'(x, y) + Iy on 
Vo has eigenvalue one. One easily checks that U'(x, y) is unitary. Since W'(0, 0) = 
the claim now follows from the case considered first. 

Let Cj(y) — U(xj(y), y) and let Pj(y) be the projection to KerCj(y). Since Cj is 
analytic in y, its eigenprojections are analytic for y ^ (but near zero) and extend 
analytically to y — (see pJj), so this is in particular true for Pj. □ 
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